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Introduction and Work Description 

This final report covers the work conducted by Brigham Young 
University personnel in the analysis and interpretation of data 
obtained by the Orbiter Retarding Potential Analyzer (ORPA) during 
the funding period of October 1, 1990 through September 30, 1994. 
During this 4 year period $46,130 was received, of which $12,000 was 
spent on faculty wages and $18,400 on student wages. The students 
working on this were Steven McNeil, Laralee Gordon, Dawn Gifford and 
Noe Yamaguchi. Steven and Laralee (1991, 1992) were graduate 
students who focused their efforts primarily on the data reduction, - 
specifically writing programs used to prepare the data (.TAB) files 
for plotting and analysis. Dawn Gifford (1993) and Noe Yamaguchi 
(199 4) were honors students who used analyses of the ORPA 
suprathermal electron data as the focus of their honors theses. 

Much of the early period of this grant was spent reanalyzing 
the ORPA data in response to the three revisions of the primary 
reduction program, LSFK, that were received approximately one year 
apart. The first revision of this program (7273 lines of code) was 
received from NASA-Ames on 11 September, 1991. The second revision 
was received 12 Aug 1992 (7 653 lines of code) and the third on 3 
September 1993 (7827 lines of code) . Between revisions 1 and 2 
there were two changes sent via e-mail in subroutine EPS (27 Nov, 
1991 and 6 Dec 1991) , and one change each to subroutines IONPK and 
SPLREA. The major difference between the 2nd and 3rd revisions was 
the extension of the analysis interval of the ORPA data to full 
orbits. We were informed in the Spring of 1993 that this expanded 
orbit analysis version was coming, but fortunately, for the specific 
program we had settled on it was not necessary to wait for revision 
3; i.e., we could still proceed with a first order determination of 
the extent to which the sprathermal electron densities derived from 
the ORPA instrument tracked the OPA proton densities in the region 
of the orbit beyond the Venus ionopause. 

During the Spring-Summer terms (typically April 20 thru August 
20) of 1993, Dawn Gifford used the second revised version of LSFK to 
determine the extent and under what conditions the ORPA and OPA 
densities tracked. Since the electrons in the solar wind are 



derived from ionization of Hydrogen, Helium and other atoms in the 
solar atmosphere, it was expected that if, beyond the ionopause, the 
suprathermal electrons were in fact solar wind electrons, then the 
two densities should agree within something approachihg 5-10%. Dawn 
found that during the times that the spacecraft potential was 
greater than 3.0 volts the two densities tracked very well. There 
is a second order effect observed at voltages above 3.0 volts, and 
with the extended orbit analysis available with the final version of 
LSFK it will be possible to determine more accurately what the form 
of the expression describing the relationship actually is 
(preliminary results suggest a simple linear relationship between 
the orbserved ORPA suprathermal electron/OPA density ratio and the 
spacecraft potential). She interrupted her studies Fall term 1993 
through Spring 1995 (LDS Mission) . Although most of the work on her 
thesis was completed prior to the date she left (September, 
1993), she just finished her thesis defense corresponding to the date 
of this report. Although dated later than that of the student who 
followed her in this project, her thesis is found at the end of 
this report as Appendix A. 

The second student, Noe Yamaguchi, used Dawn Gifford's results 
and studied 366 plots of the suprathermal electron density and 
temperature. She used the shock jump ratios to determine the 
equivalent upstream, mass-loaded mach number of the solar wind, and 
the corresponding charge exchange pickup of 0+ ions from Venus 
necesarry to explain the decrease in the mach from its defined value 
on the basis of the solar wind flow speed and magnetosonic velocity. 
A copy of her undergraduate honors thesis is included as Appendix B 

The BYU principle investigator (D. Jones) was away from the 
university from September 1994 through August 1995. As of January 1 
199 6, the work on this analysis has been started up again. A 
student is becoming familiar with the reduction and analysis 
programs and the data in anticipation of their graduate studies that 
will start this coming Fall. By that time, William Knudsen, the 
Principal Investigator of the ORPA instrument will have moved to 
Utah Valley and will be an adjunct professor in the Physics 
Department at BYU specifically helping in the analysis and 
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interpretation of the ORPA data. Hence, although the several 
modifications that occurred to the principal computer program used 
to reduce the ORPA data occurred in such a way as to essentially set 
the analysis effort back several times, as a minimum we intend to 
complete the determination of the electron density, temperature and 
velocity behind the bow shock to the point where a major paper can 
be submitted to the JGR, and,- more specifically,- to provide these 
results to theorists who need them to complete their modelling of 
the Venus bow shock and magnetosheath (see, e.g., Murawski and 
Steinolfson, 1996) . A cooperative shock analysis involving both the 
ORPA suprathermal electron and magnetometer data (Chris Russell and 
Khrishan Khurana) will also begin shortly. Hence, although the 
continual modification of the ORPA computer programs basically 
prevented us from getting the analysis to the point of being able to 
submit any papers to the JGR, nevertheless the funding has provided 
the necessary "seed money" to set an important analysis program in 
place here at BYU. 

Summary 

As a result of this analysis program it has been found that the 
ORPA suprathermal electron densities corresponding to spacecraft 
potentials of under 3.0 volts basically tracked the OPA densities 
fairly well although frequently the ORPA data appeared quite noisy. 
On the other hand, the ORPA data corresponding to spacecraft 
potentials of greater than 3.0 volts tracked very well and was 
considerably less noisy (see Figures 6,7, and 8 of Appendix A). 
Hence, using spacecraft potential dependent proportionality constant 
(near unity) , the ORPA suprathermal electron densities basically 
represent solar wind electron densities, thus allowing the high 
resolution study of Venus bow shocks using both magnetic field and 
solar wind electron data. A preliminary analysis of 3 66 bow shock 
penetrations was completed using the the solar wind electron data as 
determined from the ORPA suprathermal electron densities and 
temperatures, resulting in an estimate of the extent to which mass 
loading pickup of 0 + (UV ionized O atoms flowing out of the Venus 
atmosphere) upstream of the Venus obstacle occurred (pickup of 0 + 
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averaged 9.95%, ranging from 0.78% to 23.63%). The seed money- 
provided for this research effort will result in detailed, high 
resolution solar wind electron density and temperature data that, 
when combined with the UCLA magnetometer, will provide excellent 
data needed to constrain models presently being developed for the 
solar wind-Venus ionosphere interaction. 

Two papers have been given during the period of this research 
grant : 

D. Gifford, D. E. Jones, N. Yamaguchi, W. C. Knudsen, J. D. Mihalov, 
"Pioneer Venus Orbiter Retarding Potential Analyzer Observatins of 
the Electron Component of the Solar Wind and Observations of Venus 
Bow Shock Crossings and the Magnetosheath Region", EOS . 74., 493, 
1993. (Abstract of paper given at the Fall, 1993 Meeting of the 
American Geophysical Union, San Francisco, California, December 
1993) 

N. Yamaguchi, D. E. Jones, and W. C. Knudsen, "Pioneer Venus Orbiter 
Retarding Potential Analyzer Observations of the Electron Component 
of the Solar Wind, and of the Venus Bow Shock and Magnetosheath " , 
(paper presented at the Spring, 1994 Meeting of the Utah Academy of 
Sciences, Arts., and Letters, Ogden, Utah) 

References : 

Murawski, K. , and R.S. Steinolfson, Numerical Simulations of Mass 
Loading in the Solar Wind Interaction with Venus, JGR, 101 , 2547, 
1996. 
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Abstract 


In this study, we compared the data from two instruments on the 
Pioneer Venus Orbiter (FVO). The two data sets are the suprathermal electron 
data from the retarding potential analyzer (ORPA) and the proton data from the 
plasma analyzer (OPA). This study was designed to determine if the ORPA 
suprathermal electron data obtained outside of Venus’ ionopause represented 
solar wind electrons and could be used to increase the time resolution of the 
solar wind data from FVO. To discover whether or not this was the case, we 
compared the ORPA suprathermal electron data with the OPA solar wind 
proton data. We determined that the two sets correlate well. This indicates 
that the ORPA suprathermal electron measurements outside of the Venusian 
ionopause are in fact measurements of solar wind electrons. With this data, it 
will now be possible to conduct high time resolution studies of the bow shock 
and the region between the ionopause and the bow shock. These data will 
provide constraints on the several models currently being developed to 
represent the interaction of the mass loaded solar wind with Venus' ionosphere. 
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Introduction 


A bow shock is a paraboloid-shaped shockwave that forms around a 
planet similar to one that forms around an airplane traveling faster than the 
speed of sound. In the airplane’s case, the plane is flying at supersonic 
velocities in the air’s reference frame; the air has to move around the plane 
and forms a shockwave. In the case of a planet, the planet is traveling faster 
than the speed of sound in the solar wind’s reference frame (actually the solar 
wind is supersonic, not the planet, but the effect is the same), so the solar 
wind-planet interaction forms the bow shock. In the interplanetary medium, 
the speed of sound is typically around 30 km/ s and the solar wind speed is 
about 500 km/s. A rough diagram of a bow shock around Venus is shown in 
Figure 1. 

The bow shock can be detected by measuring several quantities, such as 
solar wind proton and electron densities, temperatures and velocities and the 
interplanetary magnetic field upstream (in front of the planet on the sunward 
side, see Fig. 1) of the planet. To get a full understanding of the bow shock, all 
of these measurements are vital. When the spacecraft taking the 
measurements passes through the bow shock, the values of the measured 
quantities will change drastically. Taking several measurements of these 
changes while passing through the bow shock at different places helps map the 
structure of the bow shock, which can change with varying conditions in the 
solar wind, planetary ionosphere, planetary magnetic field, and interplanetary 
magnetic field. The time resolution of the measurements must, however, be 
high enough to reliably measure the jumps in each quantity. 





In May 1978, NASA launched the Pioneer Venus Orbiter (FVO). FVO 
orbited Venus every twenty-four hours until Fall of 1992, taking measurements 
with a number of instruments designed to study the planet and its 
environment 1 . One of the questions FVO was devised to answer was: what is 
the structure of Venus’ bow shock and how does it respond to solar activity? 

Up to now only the magnometer could be used to study the bow shock. High 
time resolution plasma data taken while passing through the bow shock are 
needed to further understand the shock’s structure. One of the instruments on 
the orbiter, the plasma analyzer (OPA), was designed to take measurements of 
the solar wind 2 . Because the instrument was designed primarily to sample the 
upstream solar wind, the measurements closer to the planet were not taken 
often enough to adequately study the bow shock. In addition. OPA only took 
good measurements of the solar wind protons and not of the other quantities 
necessary for a complete study of the bow shock. 

Another instrument, the retarding potential analyzer (ORPA), was 
designed to take measurements of the ionosphere, including the densities, 
temperatures and velocities of ions, thermal electrons, and suprathermal. or 
high energy, electrons 2 . Because this instrument was designed to measure 


'More information about PVO is contained in "Pioneer Venus Orbiter - Ten 
years of Discovery." 

2 For a detailed description of the OPA instrument, read "The Pioneer Venus 
Orbiter Plasma Analyzer Experiment" by Intriligator. Wolfe, and Mihalov. 

3 A complete description of the ORPA instrument is contained in the papers by 
Knudsen et al. and the "Detailed instrument Description of the Pioneer-Venus 
Orbiter Retarding Potential Analyzer" prepared by Knudsen. 



near the planet the measurements were taken much more frequently than 
those taken by the OPA instrument 

Dr. William Knudsen, the principle investigator of the ORPA 
instrument has proposed that the suprathermal electrons ORPA measured 
when the orbiter was outside the Venusian ionosphere were actually solar wind 
electrons. If Dr. Knudseris idea is correct, when FVO was outside the 
Venusian ionopause. the outer boundary of the ionosphere (see Fig. 1), ORPA 
became a solar wind electron measuring device when operating in the supra- 
thermal electron mode (Knudsen Private Communication with Dr. Jones). If 
this is the case, because the OPA instrument measured only proton data and 
these only with low time resolution, the high time resolution ORPA supra- 
thermal electron measurements are critical for a better understanding of the 
region at and behind the Venusian bow shock and its dependence on solar 
activity. The main question is then: can the ORPA suprathermal electron 
measurements be used as reliable high time resolution solar wind electron 


data? 



Background 


As stated above, one of the factors PVO was designed to investigate was 
the structure of the bow shock around Venus. As different supersonic speeds 
of an airplane form different shapes and structures of shockwaves, changing 
solar wind speed, changes in solar activity, different conditions in the 
ionosphere, differences in the planetary magnetic field, and changes in the 
interplanetary magnetic field also modify the characteristics of a planetary bow 
shock. Venus’ bow shock is particularly important to study because Venus, 
unlike the Earth, has no magnetic field of its own. Earth's bow shock has been 
widely studied and as a result is relatively well understood. The lack of a 
Venusian magnetic field and differing conditions in the Venusian ionosphere 
form a very different bow shock around Venus, one that is not well understood. 

Earth’s strong magnetic field causes the Terrestrial bow shock to form 
much farther from the planet than Venus' does (typically around ten Earth 
radii from the earth compared to a little more than one Venus radius from 
Venus), keeping the solar wind well away from the Earth's ionosphere where 
particles could otherwise be picked up by the solar wind. Venus has no 
magnetic field, so the solar wind interacts with particles, particularly heavy 
ions, from the Venusian ionosphere. The solar wind picks up these heavy ions, 
causing mass loading which reduces the solar wind’s Mach number, the ratio of 
solar wind velocity to sound speed in the medium, by a factor of three to five. 
The amount of mass loading that occurs, and in turn the amount of solar wind 
slowing, depends upon solar activity and the resulting conditions in the 
Venusian ionosphere. 
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All of these factors form a very different bow shock around Venus than 
around the Earth. Study of this bow shock would greatly enhance our 
understanding of the structure of this and other bow shocks around non- 
magnetic planets and satellites. In addition, the process of bow shock 
formation around a planet occurs only under collisionless (relatively few 
collisions in the solar wind in transit from the sun to Venus) conditions, which 
are unattainable in the laboratory. Any increase in information regarding this 
process will greatly add to our understanding of the physics of collisionless 
plasmas. 

Measurements like those taken of the solar wind by OPA could 
contribute significantly to the study of Venus bow shock and its associated 
physics if they were of high enough time resolution. The highest time 
resolution of the OPA data is, however, only about eight minute intervals. 

Since the time during which the arbiter passes through the bow shock is only a 
few minutes or less, the spacecraft would pass through the bow shock before 
OPA could measure any details. 

The ORPA instrument, however, can yield several suprathermal electron 
measurements per minute while crossing through the bow shock. According to 
Dr. Knudsen's suggestion, the suprathermal electron measurements outside of 
the ionopause actually represent solar wind electrons. If so, the ORPA 
suprathermal electron data would be solar wind data of sufficient time 
resolution with which to study the bow shock. 

To test Dr. Knudsen's idea, we can compare the OPA proton densities 
with the ORPA suprathermal electron densities. In the solar wind, the electron 
and proton densities are almost identical, differing only by a small amount. In 
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the sun, where the solar wind originates, the extreme temperatures ionize all 
elements. The vast majority of atoms in the sun are hydrogen, which is only a 
proton when ionized. Helium forms the next largest fraction of the mass in the 
sun and other elements only a tiny portion. When the elements other than 
hydrogen are ionized by the heat of the sun. like hydrogen they form electrons, 
but the resulting ions are not single protons. Thus the electron density in the 
solar wind will be greater than the proton density by the amount of non- 
hydrogen ions formed in the sun (about 4%). It is simple, however, to 
compensate for this difference by accounting for the amount of helium ions 
likely to be in the solar wind. The fraction of other ions is negligible. Because 
the difference in densities is so small, the ratio of proton density to electron 
density in the solar wind is close to one. Thus, if the ratio of OPA proton 
density to ORPA suprathermal electron density is close to one. this indicates 
that Dr. Knudsen’s ideas are correct: outside of the ionopause the 
suprathermal electrons ORPA measured are indeed solar wind electrons. 

Yet another prospective benefit of this study is bow wave measurements. 
A bow wave is similar to a bow shock, but the changes in plasma parameters 
when passing through it are not as dramatic. It is formed when the mass 
loading upstream of an obstacle is slow enough and over a long enough 
distance that it does not produce the abrupt changes in solar wind speed and 
other quantities that indicate a bow shock. This occurs mostly upstream of 
comets, where the low gravity allows the neutral particles to boil off out to great 
distances upstream of the comet’s nucleus. Here they are ionized and interact 
with the solar wind, causing mass loading to occur over a very large distance. 
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On a graph, a bow shock looks like a near-vertical cliff connecting two 


very different levels of flat ground. In contrast a bow wave looks like a 
relatively gentle incline connecting two more similar heights. Figure 2 shows 
rough models of density c hang es when passing through a bow shock and a bow 
wave. The slope of the incline depends on the speed of the solar wind, which in 
turn depends on how many heavy ions it has picked up. 

Because of Venus’ unique environment, particularly the high amount of 
mass loading that occurs, the shockwave may occasionally not be abrupt 
enough to be a bow shock, but only a bow wave. Currently the locations of the 
bow shock have been determined by magnetic field data. These counted only 
abrupt changes as any form of shock: in essence, they have been used to detect 
only bow shocks, not bow waves. High time resolution data obtained from the 
ORPA instrument would make it possible to locate and study a combined total 
of between 5,000 to 10.000 bow wave and bow shock crossings. This would 
contribute to a much more accurate picture of the solar wind-Venus interaction 
and ultimately to our understanding of collisionless shock formation under 
varying conditions. 

The purpose of this study was to analyze and compare the OPA proton 
and ORPA suprathermal electron data to determine if the two sets correlate well 
enough to indicate that ORPA suprathermal electron data can be used as solar 
wind data. 
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Figure 2, Sample Density vs. Altitude Corves Passing Through Bow Shock (a) 
and Bow Wave (b). 



a. Bow Shock 



Orbiter Altitude Decreasing (inbound) 



b. Bow Wave 
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Methodology 


Analyzing the data involved a long and cumbersome process to get the 
data into useable form. Figure 3 shows a flow chart of the process. The ORPA 
data comes from NASA-Ames Research Center on large tape reels in reduced 
data records or "RDR" files for each orbit. We first loaded these files onto a 
VAX system and used a program produced by Ames to put the data into 
useable form. We then ran a program to tabulate the data into readable tables. 
Next we used a simple program to divide these files into separate files for each 
ORPA mode. i.e. thermal electron, ion. and suprathermal electron, and stripped 
the labels from the tables to make the data plottable. 

The suprathermal electron data files then needed to be rearranged into a 
format compatible with the OPA data, which we received from UCLA. First we 
modified the ORPA date/time tag to match the OPA tag. The OPA data files 
have a year-month-day-hour-minute-second date/time tag. Since the ORPA 
data comes as individual orbit data files, the files have no date tag. In addition, 
the time tag is a single number combining the hour, minute, and second. For 
example, the data for the tenth hour, twelfth minute, and 59.2 second would be 
tagged 101259.2. These incompatibilities, however, are easy to correct. We 
wrote a simple program to input the date, which we obtained from a table of 
periapsis dates for each orbit, and convert the time tag to hour-minute-second 
format for each ORPA file. 

Another difference between the OPA and ORPA data is that the OPA data 
comes in one-hour averages. In order to compare the two data sets, we 
developed another simple program to calculate hourly averages from the ORPA 




data, using the same start and stop times for each average that were used to 
obtain the OPA averages. Once the ORPA data was in one-hour averages, the 
average files, which were still in individual files for each orbit, needed to be 
appended into one long file and combined with the OPA data, synchronizing 
each ORPA average with the corresponding OPA average. 

To determine if the data sets match well enough, we made scatter plots 
of the OPA-ORPA data file, plotting OPA proton density versus ORPA supra- 
thermal electron density. If they were exactly the same, this would produce a 
straight line with a slope of one and a y-intercept of zero. This, however, would 
not occur since the densities are necessarily not identical because of the extra 
ions in the solar wind. A plot of a good match would roughly outline a straight 
line with a slope close to one and an intercept close to zero. 

As well as giving a general idea of how well the data match, the plots 
indicated what data might be inaccurate. Any data well outside of the main 
group could be questionable. We used the plots to find the data points which 
fell significantly outside the main group and then analyzed the points to 
determine why they were anomalous. 

After eliminating questionable data, we ran a linear least squares fitter 
on the data file, which fits a line to the data, giving the slope and intercept of 
this line as well as the correlation coefficient, which tells how well the data fit 
the line. The slope of the line determines the proportionality constant that 
relates the two data sets. A high correlation coefficient, a slope close to one, 
and an intercept close to zero indicate that the data match well and essentially 
that Knudsen’s ideas are correct; the ORPA suprathermal electrons are solar 
wind electrons. 
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Analysis and Results 


We processed the ORPA suprathermal electron data for orbits 1-120. A 
scatter plot of these versus the OPA proton data produced the predicted 
straight line with a slope close to one and a y-intercept relatively close to zero. 
Figure 4 shows the plot of these data and the best linear fit. The linear fit 
indicated that this line had a slope of .740, an intercept of 5.598, and a 
correlation coefficient of .631. Though the slope was relatively close to one. the 
correlation coefficient and the intercept were not good enough to represent a 
good match. 

There were, as can be seen in the figure, several data points that were 
well outside the main group. Searching through the data, we tried to determine 
the reason these points were anomalous. Any ORPA averages in which the 
variance divided by the average value was over .15 had already been 
eliminated, removing much of the data that was very far from the norm. In 
addition, any data which was taken below an altitude of 2700 km had been 
removed. This should exclude any data taken inside the ionopause (only on the 
sunward side, not in the ionotail), which would be invalid for this study. Any 
OPA averages with less than four measurements per average and any ORPA 
data where the value for the thermal electron density was zero were also 
excluded. 

Examination of the data determined that the points well outside the 
main group were probably not good data. Statistics alone indicates that any 
data far from the dominant group are probably questionable and should be 
discarded. We first tried eliminating any ORPA averages that had been 
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Figure 4. OPA Proton Density vs. ORPA Supratlicrmal Electron Density (#/cm 3 ). 

Straight line shows best linear curve fit. ORPA data taken below 2700 km 
and OPA data with less than four samples per average are excluded. 
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computed from data that were not spread uniformly throughout the hour 
interval, removing all averages in which any two measurements in the hour 
were over ten minutes apart This made certain that there were at least five 
measurements per average and that these spanned the entire hour. Figure 5 
shows the scatter plot of the improved data set. This plot looked considerably 
cleaner than the last though it still had some points away from the group. The 
linear fit produced a line with a slope of .999, a y-intercept of 1. 175, and a 
correlation coefficient of .665. These figures were much better than those from 
the last data set. There were still, however, some stray data points. 

Subsequent conversations between Dr. Jones and Dr. Knudsen and a 
closer examinations of some of the data files indicated that these questionable 
data points were related to the effect of spacecraft potential on the ORPA 
instrument. Plotting the ORPA suprathermal electron density versus the 
orbiter potential for several orbits indicated a definite trend. Figure 6 shows 
such a plot for orbit 9. Above a potential of 3, the densities are relatively 
closely spaced, but at a potential of 2.64 the densities scatter. Figures 7 & 8 
show suprathermal electron pressures throughout orbits 22 and 40. There are 
distinct lines representing different vehicle potentials. The upper trace 
corresponds to 2.64 volts, the middle to 3.74 and the lower to 5.02 and 6.48 
volts. Clearly, an algorithm can be developed to relate the spacecraft potential 
to density and eliminate the discrepancies. In any case, in the shock 
equations, it is the ratio of densities before and after the shock that matters 
and this ratio would be the same as long as we use the data corresponding to 
one vehicle potential. 
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Figure 5. OPA Proton Density vs. ORPA Supratliermal Fleet ron Density (H/cm 1 ). 

Straight line shows best linear curve fit. ORPA data with over ten minute intervals 
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Figure 6. ORPA Suprallicrmal Electron Density vs. Spacecraft Potential, Orbit 9 
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Figure 7. ORPA Supralhermal Electron Pressure vs. Altitude, Orbit 22. 



8. ORPA Suprathermal Electron Pressure vs. Altitude, Orbit 40. 



To investigate the effect of vehicle potential on the data, we plotted three 
different graphs of the data in orbits 1-79, one with all data points included 
regardless of potential, one with only those with potential greater than 2.0 volts 
included, and one with only those with po tent i a l greater than 3.0 volts 
included. The only restrictions on these data sets were: altitude greater than 
2700 km, suprathennal electron density greater than zero, number of OPA 
samples per average greater than three, and the given spacecraft potential 
limits. Those with potentials of under 3.0 volts included had a large amount of 
data points clustered at the bottom of the plot (see Fig 9 & 10). The one with 
only potentials greater than 3.0 volts looked much better (see Fig. 11). As 
noted previously, this indicated that vehicle potentials of less than 3.0 volts 
adversely affected the data. 

The plot of the data with spacecraft potential greater than 3.0 volts still 
had some problem points, however, and suggested that elimination of 
potentials less than 3.0 volts would not take care of all the orbits. More 
examination of the data showed a cyclic effect in some of the orbits where the 
ORPA and OPA data didn't match. The potential would alternate (not in a 
periodic cycle, but still visibly cyclic) between two or three different potential 
levels, and the density values would change dramatically with the potential. 
During a visit to BYU, Dr. Knud sen indicated to Dr. Jones that this could be 
because of the orientation of the ORPA sensor. When the sensor axis vector 
has any component towards the sun the solar extreme ultraviolet radiation 
(EUV) on the sensor produces extra suprathermal electrons and thus affects 
the density readings. This means any data where the sensor axis vector has 
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Figure 10. OPA Proton Density vs. ORPA Suprathermal Electron Density (#/cm 3 ). 

Straight line shows best linear curve fit. OPA data with less than four sc 
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Figure 12. OPA Proton Density vs. ORPA Suprathermal Electron Density (#/cm 3 ). 
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any component towards the sun needs to be eliminated in a further study. 

This could clean up the data well enough to get an accurate proportionality 
constant 

For the purpose of this study, however, we simply eliminated all data 
points in the most questionable orbits: 13, 14 and 25, The scatter plot of this 
data. Figure 12, looked much cleaner. The data clearly outlined a straight line 
with a slope close to one and an intercept close to zero. The linear least 
squares fit showed that this line had a slope of .842, an intercept of -.036, and 
a correlation coefficient of .756, The correlation coefficient was significantly 
higher than the previous fits and the intercept was very nearly zero. There 
were still some data points away from the main group, but they were on both 
sides of the curve fit, balancing each other out well enough to produce a 
relatively accurate preliminary proportionality constant. 

This result strongly supports Dr, Knudsen's idea that the ORPA 
suprathermal electron data outside of the ionopause does indeed represent 
solar wind electron data. The preliminary proportionality constant relating the 
OPA proton data and the ORPA suprathermal electron data is .842, the slope of 
the line. Multiplying the ORPA suprathermal electron density data by this 
number would produce high resolution solar wind electron densities with which 
to study the Venusian bow shock region. 
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Figure 11. OPA Proton Density vs. ORPA Supratliermal Electron Density (#/cm 3 ). 
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Representative Bow Shock Crossings 


To illustrate the high resolution bow shock data which ORPA gives, we 
plotted some bow shock crossings from different orbits. Figure 13 shows an 
inbound bow shock crossing during the first orbit. In this plot, the step up in 
suprathermal electron density as the vehicle passes through the bow shock is 
clearly depicted. Figure 14 shows an outbound crossing during orbit 61. The 
step is also clear in this figure, but the density decreases closer to the planet. 
Figure 15 is a plot of an outbound crossing during orbit 78. This shock is 
much more dramatic than the last two, with density increasing to about 150 
per cubic centimeter as opposed to around 40 per cubic centimeter in the last 
two. Orbit 78’s inbound crossing (Fig. 16) shows an even higher density jump, 
up to about 170 per cubic centimeter. 

Both Figures 14 and 15 show multiple density levels in the shock jumps 
caused by varying spacecraft potential. This indicates that much of the scatter 
remaining in Figure 12 is probably due to the vehicle potential. At a fixed 
potential, the suprathermal electron densities will undoubtably track the solar 
wind proton densities exceptionally well. 

Figure 17 illustrates the difference in sampling rate between the OPA 
and the ORPA data. The dotted vertical lines show approximately where the 
eight-minute OPA measurements would occur during orbit 78’s inbound bow 
shock crossing. It is clear that the OPA time resolution is much less than the 
ORPA and is not high enough to obtain an accurate picture of the jump in 
plasma parameters at the bow shock which occurs typically in a fraction of a 
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Figure 13. Bow Shock Crossing Orbit 1 Inbound 
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Figure 14. Bow Shock Crossing Orbit 61 Outbound 
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Figure 15. Bow Shock Crosang Orbit 78 Outbound 
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Figure 16. Bow Shock Crossing Orbit 78 Inbound 
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Figure 17. Bow Shock Crossing Orbit 78 Inbound 

Dotted lines show approximate OPA eight-minute intervals. 
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minute. These plots illustrate the high resolution data which the ORPA 
instrument produces and that this data is a great improvement over the OPA 
data for studying the bow shock. 
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Conclusions 


This study strongly indicates that the ORPA suprathermal electron data 
taken outside of the Venusian ionopause represents solar wind electron data. 
Our analysis showed that the ORPA suprathermal electron density data 
matches the OPA proton density data veiy well. The preliminary proportion- 
ality constant relating the two data sets is .842. The data produced using this 
constant would be high time resolution solar wind electron density data with 
which to study the Venusian bow shocks and bow waves. This is not possible 
with the low time resolution solar wind proton data currently available from the 
OPA instrument. 

A future study is needed to determine a more accurate proportionality 
constant by eliminating any densities which may have been effected by solar 
EUV, accounting for the relationship between density and vehicle potential, 
plotting more orbits (probably up to orbit 1500), and comparing the ORPA data 
with the OPA eight minute measurements, not just the hour averages. Though 
this further work is clearly needed, the present study strongly supports Dr. 
Knudsen's suggestion that when outside the ionopause. ORPA suprathermal 
electron data represents solar wind electron data. With this high time 
resolution solar wind data, potentially over 5,000 bow shock penetrations can 
now be studied under a full cycle of solar EUV variations, leading to the most 
comprehensive study of the interaction of the supersonic solar wind with a 
nonmagnetic planet that has ever been attempted and thus dramatically 
increasing our understanding of the physics of shock formation in a 
collisionless medium. 
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Appendix of Programs 


(Note: This appendix does not contain all of the 
programs used, just some of the shorter ones we developed.) 



PHOTODATA.FOR 


This program asks the user to input the year, month, and day of the orbit 
file to be processed. It then adds this date tag onto the data, converts the ORPA 
time tag to the OPA hour-minute-second format, and writes the newly formatted 
data to a new file. To run this program, the user must assign the data file to be 
processed as FOROOl and the name of the new file as FOR002. 
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DIMENSION XA(30000) ,01(30000} ,T1 (30000) ,SZA(30000) ,EP1 (30000) , 
*D2 (30000) ,T2 (30000) ,EP2 (30000) 

CHARACTER TITLE* 3 5 
83d 

888 FORMAT (IX, 'ENTER YEAR, XYR (14)') 

ACCEPT 889, IYR 

889 FORMAT ( 14 ) 

TYPE 890 

890 FORMAT (IX, 'ENTER MONTS, MON (12) ') 

ACCEPT 891, IMON 

891 FORMAT ( 12 ) 

TYPE 892 

892 FORMAT (IX, 'ENTER DAY OF MONTH, NDAY (12)' ) 

ACCEPT 893, IDAY 

893 FORMAT (12) 

IM=1 

C ELECTRON PRESSURE WILL BE IN UNITS OF DYNES/ CM2 TIME E10. 

ALPHA=1 . 6E-02 

11 READ (1,333, END- 77.3 4 ) N1 , P2 , P3 , P4 , P5 , PS, P7 , P8 , P9 , P10 , Pll , P12 , P13 
333 FORMAT (15 , 2F9 .1 , F3.1, F10 .2 , 2F11 . 3 , SF10 . 3) 

IF (P6 .LE. 0.0) GO TO 11 
IF (P8 .LE. 0.0) GO TO 11 
IF (P10 .LE. 0.0) GO TO 11 
IM-IM+1 
TT=?2 

IHR=INT (TT/10000 . ) 

DELT=TT- 10000 . *REAL (IHR) 

IMIN-INT (DELT/100 . 0 ) 

SEC=TT- (10000 . *REAL (IHR) +100 . *RSAL (IMIN) ) 

XA(IM) =P3 
SZA(IM) *?4 
D1 ( IM) =?S 
Tl { IM) =?8 
D2 (IM) =P10 
T2 ( IM) =P12 

EP1 (IM) =ALPHA*D1 ( IM) *T1 ( IM) 

EP2 (IM) =ALPKA*D2 (IM) *T2 (IM) 

WRITE (2 ,3434) IYR, IMON, IDAY, IHR, IMIN, SEC, XA(IM) ,SZA(IM) , 

*D1 (IM) , Tl ( IM) , S?1 (IM) ,D2 (IM) ,T2 (IM) ,EP2.{IM) 

3434 FORMAT (15 , 413 , F5 . 1 , IX, F8 . 1, FS .1, F11.2 , 5F11. 3) 

GO TO 11 
7734 CONTINUE 

CLOSE { UNIT- 1) 

12 CALL EXIT 
END 
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PHOTOAVE 1 J?OR 


This program takes one-hour averages of the data file assigned FOROO 1 and 
writes them to the file assigned F0RQ02. It puts high junk numbers in as the 
data for any averages with less than three measurements per average and removes 
any data taken below 2700 km. 



noon noon 


THIS PROGRAM REQUIRES PHOTODATA. DAT ASSIGNED AS FORGO 1 AND 
PHOTOAVE.DAT ASSIGNED AS FOR002 

DIMENSION DENS (9000) , TEMP {9000} ,NT(9000) 

TYPE 80 

80 FORMAT (IX, 'TAILBX POLARITY : P - BX > 0,N - BX < 0 , A® ALL ' ) 
ACCEPT 81, SIGN 

81 FORMAT (Al) 

1003 CONTINUE 
KX-1 

1002 READ (1,3434, END-900) IYR, IMQN, IDAY , IHR, BUN, SEC , XA, SZA, 

*D1 , T1 , EP1 , D2 , T2 , EP2 

3434 FORMAT ( 15 , 413 , F5 .1 , IX, F8 . 1 , F6 . 1 , Fll . 2 , 5F11 . 3 ) 

IF (XA .LT. 2700.0) GO TO 1002 
NT (IOC) -IHR 

IF (IHR -GT. NT(1) ) GO TO 1004 
IF ( IHR .LT. NT(1) ) GO TO 1004 
DENS (KX) -REAL (D2) 

TEMP (KK) -REAL (T2) 

KX-KX+1 
10 GO TO 1002 

1004 CONTINUE 
KPTS-KX-1 

IF ( X?TS .LE. 2) GO TO 18 

CALL XFIT (DENS, KPTS, 0, DENAVE, SIGDN, SIGD) 

CALL XFIT (TEMP , KPTS , 0 , TEMAVE , SIGTN, SIGT) 

GO TO 19 

18 DENAVE-0.0 
SIGDN-999 . 9 
SIGD-999 . 9 
TEMAVE =0 - 0 
SIGTN-999 .9 
SIGT-999.9 

19 CONTINUE 

WRITE (2 ,3800) IYR, IMON, IDAY, NT (1) , KPTS, DENAVE, SIGDN, SIGD, 
♦TEMAVE , SIGTN, SIGT 
3800 FORMAT (IX, I4,4I3,6F10.2) 

SDENS-REAL (D2 ) 

STEMP-REAL (T2 ) 

NT (1) =NT(KK) 

KX-I 

DENS (KX) -SDENS 
TEMP ! KX) -STEM? 

KX-KK-1 
GO TO 1002 
900 CONTINUE 
KPTS-KK-1 

IF (KPTS .LE. 2) GO TO 28 

CALL XFIT (DENS , KPTS , 0 , DENAVE , SIGDN, SIGD) 

CALL XFIT (TEMP , KPTS , 0 , TEMAVE , SIGTN, SIGT) 

GO TO 29 

28 DENAVE-0.0 
SIGDN-999. 9 
SIGD-999. 9 
TEMAVE-0 . 0 
SIGTN-999 .9 
SIGT-999.9 

29 CONTINUE 

WRITE (2,3800) IYR , IMON , IDAY , NT { 1 ) , KPTS , DENAVE, SIGDN, SIGD , 
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♦TEMAVE, SIGTN, SIGT 
CLOSE (UNIT»1) 

CALL EXIT 
END 


q******* *************** ******************************** 

c 

********************************* 

SUBROUTINE XFIT {X , NPTS , MODE , XMEAN, SIOfflM, SIGMA) 
DOUBLE PRECISION SUM, SUMX, WEIGHT, FREE 
DIMENSION X(l) 

C 

C ACCUMULATE WEIGHTED SUMS 

C 

IF (NPTS .LE. 2) GO TO 888 
11 SUM = 0, 

SUMX « 0, 

SIGMA = 0. 

SIGMAM = 0. 

20 DO 32 I = 1, NPTS 

21 IF (MODE) 22, 22, 24 

22 WEIGHT = 1. 

GO TO 31 

24 WEIGHT = 1. 

31 SUM = SUM + WEIGHT 

32 SUMX = SUMX + WEIGET*X(I) 

C 

C EVALUATE MEAN AND STANDARD DEVIATIONS 

C 

41 XMEAN = SUMX/SUM 

51 DO 52 1=1, NPTS 

52 SIGMA = SIGMA + (X Cl) -XMEAN) **2 
FREE = NPTS - 1 

54 SIGMA = DSQRT ( SIGMA/ FREE) 

61 IF (MODE) 62, 64, 66 

62 SIGMAM = DSQRT( XMEAN/ SUM) 

GO TO 70 

64 SIGMAM = SIGMA / DSQRT (SUM) 

GO TO 70 

66 SIGMA = DSQRT ! 1 . /SUM) 

GO TO 889 

888 WRITE (4, 9033) 

9033 FORMAT (IX, 'NOT ENOUGH DATA') 

889 CONTINUE 
70 RETURN 

END 
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PHOTOAVE2JOR 


This program is the same as Photoavel.for, but it excludes any averages 


with time jumps of over ten minutes between any two measurements. 



onnn 


THIS PROGRAM REQUIRES PHOTODATA. DAT ASSIGNED AS FOROOl AND 
PHOTOAVE.DAT ASSIGNED AS FOR002 

DIMENSION DENS (9000) , TEMP (9000) , NT (9000) , MDIFF (9000) 

C TYPE 80 

C 80 FORMAT (IX, 'TAILBX POLARITY : P « BX > 0,N . BX < 0 , A* ALL ' ) 
C ACCEPT 81, SIGN 

C 81 FORMAT (Al) 

1003 CONTINUE 
KX-1 
KD-1 
MD-0 

1002 READ (1,3434, END-900) IYR, IMDN, ID AY, IHR, IMIN, SEC, XA, SZA, 

*D1 , T1 , EPl , D2 , T2 , EP2 

3434 FORMAT {15 , 413 , F5 .1, IX, F8 .1, F6.1, F11.2, 5F11.3) 

IF (XA .LT. 2700.0) GO TO 1002 
NT (XX) -IHR 

IF (IHR .GT. NT(1) ) GO TO 1004 
IF (IHR .LT. NT(1) ) GO TO 1004 
MDIFF (XD) -IMIN-MD 
MD— IMIN 
XD-XD+1 

DENS ( XX) -REAL (D2) 

TEMP (XX) -REAL (T2) 

KX-XX+1 
10 GO TO 1002 

1004 CONTINUE 
MDIFF (XD) -59-MD 
KPTS-XX- 1 
MCOUNT = XD - XPTS 

111 IF (MDIFF (MCOUNT) .GT. 10) GO TO 222 
MCOUNT-MCOUNT+1 

IF (MCOUNT .GT. XD) GO TO 112 
GO TO 111 

112 IF (XPTS .LE. 2) GO TO 18 

CALL XFIT ( DENS , XPTS , 0 , DENAVE , SIGDN , SIGD ) 

CALL XF IT ( TEMP , XPTS , 0 , TEMAVE , S IGTN , S IGT) 

GO TO 19 

18 DENAVE-C . C 
S IGDN-59S.S 
SIGD-999 .9 
TEMAVE =C .0 
S IGTN 5*999 . 9 
SIGT-999.9 

19 CONTINUE 

WRITE (2,3 800) IYR, IMON, IDAY,NT(1) , XPTS, DENAVE, SIGDN, SIGD, 

* TEMAVE , S IGTN , S IGT 
3800 FORMAT (IX, 14, 413, 6F10. 2) 

222 SDENS-RSAL (D2) 

STEMP-REAL (T2 ) 

NT (1) -NT (XX) 

XX- 1 

DENS (XX) -SDENS 
TEMP (XX) -STEMP 
KX-XX+1 
KD— KD+1 
MD-0 

MDIFF (KD) -IMIN-MD 
MD-IMIN 
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KD-KD+1 
GO TO 1002 
900 CONTINUE 

MDIFF (KD) *59 -MD 

KPTS-KK-1 

MCOUNT»KD-KPTS 

211 IF (MDIFF (MCOUNT) .GT. 10) GO TO 322 
MCOUNT-MCOUNT+1 

IF {MCOUNT .GT. KD) GO TO 212 
GO TO 211 

212 IF (KPTS .LE. 2) GO TO 28 

CALL XFIT ( DENS , KPTS , 0 , DENAVE , SIGDN , SIQ3) 

CALL XFIT (TEMP, KPTS, 0 , TEMAVE, SIGTN, SIGT) 

GO TO 29 

28 DENAVE-0.0 
SIGDN-999.9 
SIGD-999.9 
TEMAVE-0 . 0 
SIGTN-999.9 
SIGT-999.9 

29 CONTINUE 

WRITE (2, 3800) IYR, IMON, IDAY,NT (1) , KPTS, DENAVE, SIGDN, SIGD, 
♦TEMAVE, SIGTN, SIGT 
322 CLOSE ( UNIT- 1) 

CALL EXIT 
END 


Q** ********************** ***** ************** *********** 

C 

********************************* 

SUBROUTINE XFIT (X, NPTS ,MODE, XMEAN, SIGMAM, SIGMA) 
DOUBLE PRECISION SUM, SUMX, WEIGHT, FREE 
DIMENSION X(l) 

ACCUMULATE WEIGHTED SUMS 

IF (NPTS .LE. 2) GO TO 888 
11 SUM = C . 

SUMX = 0. 

SIGMA * 0. 

SIGMAM = 0. 

20 DC 32 I » 1, NPTS 

21 IF (MODE) 22, 22, 24 

22 WEIGHT = 1, 

GO TO 31 

24 WEIGHT - 1. 

31 SUM - SUM + WEIGHT 

32 SUMX = SUMX + WEIGHT*X(I) 

EVALUATE MEAN AND STANDARD DEVIATIONS 

41 XMEAN = SUMX/SUM 

51 DO 52 1=1, NPTS 

52 SIGMA = SIGMA + (X (I) -XMEAN) **2 
FREE = NPTS - 1 

54 SIGMA = DSQRTt SIGMA/ FREE) 

61 IF (MODE) 62, 64, 66 

62 SIGMAM = DSQRT (XMEAN/ SUM) 
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GO TO 70 

64 SIGMftM - SIGMA / DSQRT (SOM) 
GO TO 70 

66 SIGMA - DSQRT (1. /SUM) 

GO TO 889 

888 WRITB{4, 9033} 

9033 FORMAT (12, 'HOT ENOUGH DATA') 

889 CONTINUE 
70 RETURN 

END 



ORPAOPAJTOR 


This program correlates the OPA proton data and the ORPA suprathennal 
electron one-hour averages, matching the data by date/time tags. -It also removes 
any ORPA data where the variance divided by the average value is greater than 
.15. It can also remove any data where the OPA density is over 60 or the ORPA 
temperature average is over 250000, but these lines are commented out in this 
version. The program requires the OPA data file assigned as FOROOl. the ORPA 
averages data file assigned as FOR002, and the output file name assigned as 
FOR003. 
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REAL NREC 

1002 OPEN (UNIT-2 , STATUS- 'OLD' , FORM- ' FORMATTED ' , RECL-11 , 

♦RECORDTYPE- ' FIXED' , IOSTAT-ISTAT, READONLY) 

12 READ (1,302, END-900) NYR, NMON, NDAY, NHR, NMIN, SEC, NREC, TEMP, 

♦DENS, VEL.AZ, POL 

302 FORMAT(I5,4i3,F7.3,F8,0,F12.Q,F9.3,F7.1,2F9.3* 

IF (NREC .LT. 5) GO TO 12 

11 READ (2, 3800, END-800) IYR, IMQN, IDAY, IHR, KPTS , DENAVE, SIGDN, SIGD , 
♦TEMAVE , SIGTN, SIGT 
3800 FORMAT (14, 413, 6F10. 2 5 
C IF (DENS ;GT. 60.0) GO TO 800 
. TEMAVE-TEMAVE*11594 .2 

IF (DENAVE .EQ. 0.00) GO TO 11 
IF (SIGDN/DENAVE .GT. 0.15) GO TO 11 
C IF (TEMAVE .GT. 250000.0) GO TO 11 

IF ( IYR . EQ . NYR . AND . IMON . EQ . NMON . AND . IDAY . EQ . NDAY . AND . IHR . EQ . NHR ) 
♦GO TO 10 

IF ( IYR . NE . NYR . OR . IMON . NE . NMON . OR . IDAY . NE . NDAY . OR . IHR . NE . NHR ) 

♦GO TO 11 

10 WRITE (3 , 3010) NYR, NMON, NDAY, NHR, DENS, DENAVE, TEMP, TEMAVE 
3010 FORMAT (15, 313, 4F10. 2) 

800 CONTINUE 

CLOSE (UNIT-2) 

GO TO 1002 
900 CONTINUE 

CLOSE ( UNIT- 1) 

CALL EXIT 
END 
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DATAFDLFOR 


This program remows airy data from the correlated ORPA/OPA data file in 
which the OPA density divided by the ORPA density is greater than 1.70 or less 
than .30. This is the program I used to eliminate any questionable data far from 
the main group. It requires the ORPA/OPA data file assigned as FOROOl and the 
output file name assigned as FOR002. 
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12 READ ( 1 , 111 , END-77) NYR , NMON, NDAY, NHR, DENS. DENAVE, TEMP , TEMAVE 
111 FORMAT (15 , 313, 4F10 .2) 

IF (DENS/DENAVE .LT. 0.30) GO TO 12 
IF (DENS/DENAVE .GT. 1.70) GO TO 12 

WRITE (2,222) NYR, NMON, NDAY, NHR, DENS , DENAVE, TEMP, TEMAVE 
222 FORMAT (15, 313, 4F10. 2) 

GO TO 12 
77 CLOSE (UNIT-1) 

CALL EXIT 
END 
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UNFTTFVO.FOR 


This Is the linear least squares fitter we used to fit the ORPA density vs. 
OPA density lines. It requires the ORPA/OPA file assigned FGR001 and the 
output file name assigned FOR002. It writes the number of data points, the y 
intercept, the variance on the intercept, the slope, the variance on the slope, the 
correlation coefficient, the average ORPA density, and the average OPA density to 
the output file. 
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C THIS IS THE PROGRAM "LINFITFVO.FOR* 

C 

DIMENSION AA(lOQOO) ,BB (10000) 

NN-0 

1 READ (1, 100, END-9 00) OPADENS / ORPADENS 
100 FORMAT (14X,2F10. 2) 

NN-NN+1 

AA(NN) -ORPADENS 
BB (NN) -OPADENS 
GO TO 1 

900 CLOSE ( UNIT- 1) 

KPTS-NN 

IF (KPTS .EQ. 0) GO TO 302 

CALL LINFIT (AA, BB , KPTS , YNT,SGYNT, YSLP , SGYSLP ) YY, AVA, AVB ) 
WRITE (2, 300) KPTS , YFT , SGYNT , YSLP , SGYSLP , YY , AVA, AVB 
300 FORMAT ( 17 , 7F9 , 3 ) 

302 CALL EXIT 
END 
C 

Q**ir± **★★*★* ********** ***^ ******* ************************* *.*%■**■*** 

c 

SUBROUTINE LINFIT (X , Y, NPTS , A, SIGMAA, B , SIGMAB , R, AVX , AVY) 
DOUBLE PRECISION SUM, SUMX, SUMY, SUMX2 , SUMXY, SUMY2 
DOUBLE PRECISION XI, YI, WEIGHT, DELTA, VARNCE 
DIMENSION X(I) ,Y{1) 


ACCUMULATE WEIGHTED SUMS 


IF (NPTS .LE. 2) GO TO 500 
11 SUM=0. 

SUMX-C , 

SUMY-0 . 

SUMX2-0. 

SUMXY- 0 . 

SUMY2-C . 

21 DO 50 1=1, NPTS 
XI = X .1) 

YI = Y(l! 

WEIGHT = 1. 

41 SUM - SUM - WEIGHT 

SUMX = SUMX - WEIGHT V XI 
SUMY = SUM:’ - WEIGH7*YI 
SUMX2- SUMX2- WEIGHT*XI*XI 
SUMXY- SUMXY- WEIGHT*XI*YI 
SUMY2- SUMY2- WEIGHT *YI*Yl 
50 CONTINUE 
C 

C CALCULATE COEFFICIEINTS AND STANDARD DEVIATIONS 
C 

AVX-SUMX/NPTS 
AVY -SUMY /NPTS 

.51 DELTA = SUM*SUMX2 - 3UMX*SUMX 
IF (DELTA .EQ, 0.) GO TO 500 
A - (SUMX2*SUMY - SUMX* SUMXY) / DELTA 
53 B - (SUMXY* SUM - SUMX*SUMY } / DELTA 
C - REAL (NPTS -2) 

VARNCE = (SUMY2+ A*A*SUM + B*B*SUMX2 
* -2 . * (A*SUMY + B* SUMXY - A*B*SUMX) ) / C 
67 SIGMAA = DSQRT (VARNCE*SUMX2 / DELTA) 
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S3 SXGMAB - DSQRT (VARNCE*SUM / DELTA) 

71 E - (SUM*SUMXY - SUMX*SUMY) / 

*DSQRT (DELTA* (SUM*SUMX2 - SUMX^SUMY) ) 
500 CONTINUE 
RETURN 
END 
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PHOTORUN.COM 


This is a sample of a command file which can be used to process the data. 
It takes the data from the plottable data files to the correlated GRPA/OPA data 
files. To use it simply change the file names and the date entered to fit the 
orbit(s) to be processed. 
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$ ASSIGN PEL0001.DAT FOROOl 
$ ASSIGN PHOTODATA01.DAT FOR002 
$RDN PHOTODATA 
1978 
12 

4 

$ ASSIGN PHOTODATA01.DAT FOROOl 
^ASSIGN PHOTOAVE 0 1 . DAT FOR002 
$RUN PHOTOAVE 

$ ASSIGN PEL0002.DAT FOROOl 
SASSIGN PHOTODATA02 .DAT FOR002 
$RUN PHOTODATA 
1973 
12 

5 

$ ASSIGN PHOTODATA02.DAT FOROOl 
$ ASSIGN PHOTOAVE 02 .DAT FOR002 
SRUN PHOTOAVE 

$ APPEND PHOTOAVE 0 2 . DAT PHOTOAVE01.DAT 
$ RENAME PHOTOAVE01.DAT PHOTOAVE . DAT 
$ ASSIGN OPA-DAT FOROOl 
$ ASSIGN PHOTOAVE . DAT FOR002 
$ ASSIGN ORPAOPA.DAT FOR003 
SRUN ORPAOPA 
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1. Abstract 


We have used the supra thermal electron data obtained by the Pioneer Venus Orbiter 
ORPA (Retarding Potential Analyzer) instrument to study the bow shock of Venus. The supra 
thermal electron density component of ORPA data was found by Gifford (1993) to track the solar 
wind proton density under certain spacecraft potential conditions. Using the Spreiter (1966) gas 
dynamic equations, we have estimated the Mach number that can be calculated from the change 
in plasma parameters that results from the interaction of the solar wind with the upper ionosphere 
of Venus. Using the Mach number calculated from Spreiter's equations and employing the mass 
loading equations developed for the solar wind interaction with comets, we have been able to 
estimate the fractional number of solar wind protons that have been replaced by atomic oxygen 
ions through charge exchange. 

The Mach number calculated had its average of 2.71 and ranged from 1.14 to 4.27 
(excluding two abnormal events). The fraction of the picked up solar wind protons that were 
replaced by 16 m p ions (oxygen ions) had an average of 9.95% and ranged from 0.78% to 
23.63%. This result takes only charge exchange and the resulting mass loading into account, 
neglecting the effect of interplanetary magnetic field pile up in the day-side ionosphere of Venus. 
At times it is clear that such is needed in order to more accurately estimate the amount of charge 
exchange pick up of 0+. 
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2. Introduction 


2.1 PVO Mission 

The Pioneer Venus Orbiter (PVO) was launched by NASA on May 20, 1978, and injected 
into Venus orbit on December 4, 1978. PVO orbited Venus every twenty-four hours until the 
Fall of 1992, taking measurements with a number of instruments designed to study the planet and 
its environment (See Figure 2.1.1). Two of the instruments on PVO designed to study plasma 
properties around Venus were the plasma analyzer (OP A) and the retarding potential analyzer 
(ORPA). These two instruments were a part of the complement of six in-situ 1 experiments 
during the mission. The data were made available to us by NASA-Ames Research Center and 
fills a vital role in the study of the plasma environment of Venus. 

Before the arrival of the PVO at Venus on December 4, 1978, no spacecraft had entered 
that planet’s ionosphere to study it as described in Tatrallyay and his coworkers’ report (1983): 

"The first spacecraft to visit Venus were either flyby, bus, or entry probe vehicles 
that provided only one or two samples of the bow shock, hence limiting the ability to 
perform statistical investigations of the physical conditions around the earth's nearest 
neighbor. ... in December 1978, the Pioneer Venus orbiter (PVO) was placed into a 24- 
hour Venus orbit and has been transmitting data almost continually since that time. One 
of the first studies to be made with the resulting data concerned the location of the Venus 
bow shock. While the identification of most crossings of the bow shock is straight 

1 “In-situ experiments measure quantities at the spacecraft's location, either by ingesting 
particles into the instruments or by measuring local, non-propagating fields, as opposed to the 
remote sensing experiments where they measure electromagnetic fields (optical, infrared, 
ultraviolet. X-ray, and radio) propagating from a reflecting or emitting source on the planet or 
within its atmosphere.*’ ( Ten Years of Discovery ) 
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Figure 2.1.1 Diagram of the Pioneer Venus Orbiter spacecraft identifying the main 
components. 



Source: Colin (1989) 
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forward in almost any of the plasma and field measurements, the subject of die location 
of the bow shock has been one of the most controversial topics of Venus studies because, 
in part, the location of the bow shock is sensitive to the physics of the interaction, and the 
physics of that interaction is at best poorly understood-" 

Numerous studies have been completed since the time PVO started receiving data on the near 
environment of Venus beginning near the end of 1978. Some of these studies involve the 
vigorous modeling of various characteristics of the bow shock. Our research follows and, at the 
same time, makes an attempt to expand these earlier investigators’ works on the Venus’ bow 
shock formation. We utilize the electron component of the PVO ORPA data, which has not been 
reported up to this point. 

2.2 Bow Shock Formation 

A bow shock is a paraboloid-shaped shockwave that forms around the planet as the solar 
wind interacts with the planet (See Figure 2.2.1). The process of planetary bow shock formation 
is similar to that which forms when a supersonic airplane passes through the upper atmosphere 
at speeds greater than that of sound. This interaction between the aircraft and the air produces a 
shock that is caused by the phenomena called steepening. Steepening occurs because the passage 
of either a fast or slow wave alters the previously undisturbed air condition and thus allows the 
trailing wave to travel faster and catch up with the first wave. 

The shock formation from the aircraft-air interaction faces much more simplified physical 
phenomena compared to what takes place when the bow shock forms from the solar wind-planet 
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Figure 2.2.1 Bow Shock Formation — Sketch of the principal regions of the terrestrial 
magnetosphere and its connections to the solar wind magnetic field 
(not in scale). 



Source: Lui (1987) 
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interaction. The latter case must deal with magneto hydrodymanic (MHD 2 ) effects which include 
die interplanetary magnetic field influences, and the complex plasma behavior of the solar wind. 
However, the basic shock formation process from the solar wind-planet interaction still shares 
the similar effects also feced by the aircraft-air interaction. When the solar wind traveling 
through space reaches the planet, it compresses both the planet’s plasma and the planet’s 
magnetic field. As the various waves of solar wind arrive at the planet, they appear to have 
similar effects as the steepening, if the amplitude of those waves are sufficient. A bow shock 
wave may be described as a thin discontinuity where the normal solar wind wave conditions are 
altered causing irreversible changes in the plasma. Such a discontinuance occurs at the bow 
shock where the plasma parameters show abrupt changes. 

We can locate the position of the interplanetary bow shock and study its characteristics by 
measuring the change of solar wind plasma 3 parameters, such as density, temperature and 
pressure. Shifts in these parameters take place during any kind of shock formations. For 
example, the bow shock formation from the aircraft-air interaction is detected by the human ear 
drum since the abrupt pressure change before and after the bow shock is large enough for the ear 
drum to sense. However, for the MHD dependent parameters, such as magnetic field, the 
discontinuance is observed only in the case of bow shock formation within the interplanetary 
medium. 

2 MHD requires use of physical laws, such as Maxwell equations, Ohm’s law, and the 
conservation of mass and momentum, to describe the characteristics of the fluid plasma -- 
Plasmas in the solar system often find themselves in unstable situations in which the MHD 
equations would predict a rapid change of configuration. 

3 Plasma is gas of charged particles in which the potential energy of a particle due to its 
nearest neighbor is much smaller than its kinetic energy 
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3. Background 


3.1 Venus- Earth Comparison 

Venus has been described as the earth's twin because of its numerous similar 
characteristics to the earth. Some of the characteristics include its mass ( 0.8 1 earth masses), 
radius (0.95 earth radii), mean density (95% that of earth), and gravity (90% that of earth). 
However, in terms of atmospheric structure and chemical composition, “Venus is anything but a 
twin of earth” (Levine 1 992). Some of the notable atmospheric structural differences are: the 
mean planetary surface temperature of Venus being 750°K, compared to about 300 3 K for earth; 
and the surface pressure on Venus being about 90 atm, while it is 1 atm for earth. Also, the most 
apparent atmospheric chemical composition difference is that “. . .carbon dioxide at 96 % by 
volume, is the overwhelming constituent in the atmosphere of Venus, while it is only a trace 
constituent in the earth’s atmosphere (0.034% by volume)” (Levine 1992). 

The most important difference between the two planets for our study lies in the 
characteristics of their planetary magnetic fields. Earth’s magnetic field is strong enough to 
stand off the solar wind pressure at about 10 earth radii. Since there is no intervention between 
the solar wind and the terrestrial atmosphere, the interaction is very “clean.” On the other 
hand, there is zero or undetectably small planetary magnetic field for the case of Venus, so that 
the solar wind can interact directly with the Venus atmosphere and ionosphere (Tatrallvay 1983). 
This contrasting magnetic nature between these two planets, which are similar in many other 
ways, urges us to perform a comparison study to further our understanding of the physics of 
plasma behavior. 
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3.2 Solar Wind Interaction with Venus 


As discussed previously, the bow shock formation in the interplanetary medium differs 
greatly from the case in the earth's atmosphere. The two cases share the similarity only in the 
context of the pressure jump occurring due to the density increase and temperature increase 
(p = nkT) at the obstacle interaction point The formation of the bow shock of Venus involves 
considerably greater complexities, resulting from the fact that the solar wind is collisionless 4 and 
is permeated by a magnetic field. 

The solar wind is a fully ionized plasma that flows continuously outward from the solar 
corona, dragging the “frozen- in” remnant of the solar magnetic field through the solar system. 
Solar wind plasma is composed almost entirely of electrons and protons, with small numbers 
(few percent) of He“. As a result of the strong outward pressure in the solar corona, the solar 
wind becomes supersonic at a few solar radii (1 solar radius = 6.97 x 10 5 km) above the visible 
surface of the Sun (the photosphere). It attains speeds in the range 250 - 750 km/s in 
interplanetary space and is believed to remain supersonic out to a distance from the Sun of 50- 
100 astronomical units (AU). Although the solar wind eventually loses its high speed as it 
encounters the interstellar gas and magnetic field, the relatively collisionless space in which they 
travel through (only 2-3 collisions between sun and earth orbit) allows them to retain their high 
speed for a long time (Levine 1992). 

Due to the absence of a planetary magnetic field in the space around Venus, the solar 

4 Although interplanetary space is generally perceived as a vacuum, it contains continuously 
changing low density plasma. Also, the solar wind is so tenuous that collisions between the 
charged particles of the plasma occur only several times on average during the entire travel from 
the sun to earth orbit (distance of 1.5 x 10 s km). Hence, the plasma is termed collisionless. 
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wind is able to interact directly with the upper atmosphere of the planet. When the supersonic 
solar wind interacts with the upper atmosphere of Venus, it produces a number of complex 
surfaces in this region such as a bow shock, an ionosheath, and an ionopause (Figure 32.1). As 
the solar wind approaches the Venus ionosphere, neutral particles from the planet move upstream 
and become subject to the interaction with the solar wind ionized particles. This interaction 
between solar wind plasma and outflowing neutral oxygen atoms from Venus results in the 
conversion of the latter into oxygen ions through charge exchange, which are then picked up by 
the solar wind. (The pick up results from any events that can cause the ionization of the particles 
in the solar wind.) Charge exchange and ionization by ultraviolet radiation from the sun are the 
two major sources of the ionization and subsequent pick up of O' during the solar wind- Venus 
interaction. 

Specifically, as the high speed hydrogen ions in the solar wind approach the ionosphere 
of the planet, they encounter the neutral oxygen atoms from the planet that are moving upstream. 
As the hydrogen ion and the neutral atoms collide with each other, charge exchange occurs as 
follows: 

H' (hot) + 0 (cold) ==> H (hot) + O' (cold), 
i.e., 

H' (150,000 K°) + 0 (300 K 3 ) 

(500 km/s) (2 km/s) 

=> H (1 50,000 K 3 ) + O' (300 K 3 ) 

(500 km/s) (2 km/s). 

Each O* formed takes the place of an H' , or a proton, of the solar wind and gets picked up 
(dragged along) by the motional electric field (E = -v x B) produced by the moving solar wind. 
This reaction is significant because the neutral oxygen atom that is flowing away from Venus at 
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Source: Colin (1989) 
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only several kilometers per second is brought up to several hundred times the original speed 
when ionized. Because the ionized oxygen weighs sixteen times more than the ionized 
hydrogen, conservation of momentum causes the slow-down of the solar wind. 

As noted above, another ionization source is the solar ultraviolet radiation. If the neutral 
oxygen atom is in the solar UV long enough, it becomes ionized by means of the reaction: 

0 + hv ==> 0 + + e\ 

Both the ionized oxygen and the electron are picked up by the solar wind flow since they are now 
subject to the magnetic field induced by the solar wind. The slow down process of the solar 
wind due to the UV radiation is the same as with the charge exchange effect. However, it is 
considered much less significant than that of the charge exchange because the ionization time by 
this mechanism is too long for the Venus case. 

The solar wind slows as it approaches the upper ionosphere of the planet. This slow- 
down of the solar wind is caused by the phenomenon called mass loading. The concept of mass 
loading is similar to that of the traffic congestion. As more vehicles come into the highway, they 
all have to slow down. Unaware of the solar wind slow down taking place at the upper 
ionosphere of the planet, the in-coming solar wind continues to send its plasma into the area, 
creating the congestion of particles at the interaction point. If the rate of mass loading of the 
solar wind is great enough, it can lead to the formation of a bow shock. This forefront section of 
the mass loading field is the bow shock. At the bow shock, the plasma parameters of the solar 
wind change abruptly. The solar wind velocity drops behind the bow shock because a large part 
of the plasma gets thermalized. As a result, the 400 km/sec purely rectilinear velocity (outward 
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from sun) becomes randomized. This shows up as a much higher plasma temperature behind the 
shock (of the order of a million degrees or so). Since die particle congestion takes place, the 
density of solar wind should also be greater. The increase in both the density and the 
temperature in the formation of the shock occur under adiabatic conditions. Consequently, more 
active interactions among particles arise, thus resulting also in a higher temperature, and 
therefore a higher pressure. 
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4. Methodology and Analysis 


4.1 Observation of the Bow Shock from the ORPA Data 

4.1.1 ORPA-OPA Data Correlation 

The measuring instrument on PVO designed primarily to study the upstream solar wind, 
the Orbiter Plasma Analyzer (OP A), was sampled so infrequently that the measurements closer 
to the planet were not taken often enough to adequately study the bow shock. In addition, the 
OPA solar wind ion temperature and electron density measurements, which would be necessary 
for a complete study of the bow shock, are not reliable. Another instrument, the retarding 
potential analyzer (ORPA), was designed to take measurements of the ionosphere, including the 
densities and temperatures of ions, thermal electrons and supra thermal (high energy) electrons. 
Because this instrument was designed to measure near the planet, the measurements were taken 
much more frequently than those taken by the OPA instrument. The time resolution of the 
ORPA data is high and the electron data is reliable. However, it was not clear how the ORPA 
supra thermal electron data related to solar wind electrons. 

In the fall of 1993, another BYU undergraduate physics student, Dawn Gifford, 
completed a study which indicates that the tracking of the ORPA and OPA density data is quite 
good most of the time if certain simple spacecraft voltage constraints are satisfied as shown in 
Figure 4.1. 1.1. As a result, because of the higher sampling rate of the ORPA compared to that of 
the OPA (comparison shown on Figure 4.1.1 .2), the former will allow a reasonably precise study 
of the location and nature of the bow shock that forms upstream of Venus in terms of the 
upstream and shocked solar wind electrons. Hence, upon carefully examining her study 
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OPA Proton Density vs. ORPA Supratberma! Electron Density (ff/cm 1 ). 

Straight line shows best linear curve fit. OPA data with less than four samples per 
average and data taken below 2700 km., ORPA data with potential less than Hire 
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of the OPA-ORPA data compatibility, we find it feasible to employ the ORPA data for our 
analysis of the bow shock. We will utilize the data available to us from the PVO mission. The 
computing programs, developed by the principle investigator of the OPA and ORPA instrument, 
have been employed to process the raw information from the PVO into supra thermal electron 
densities and temperatures. Additional computer programs have been developed to utilize only 
the supra thermal electron data for which the spacecraft potential constraints determined by 
Gifford (1993) have been satisfied. 

4.1.2 Plot ORPA Plasma Parameter Data 

Approximately 825 orbits of the ORPA density, temperature and pressure data have been 
plotted against the altitude of the spacecraft in each of the PVO orbits of Venus. By examining 
these plots, we have been able to closely study the nature of the Venus bow shock using the 
electron component of the solar wind as measured by the ORP A instrument. 

The Figures 4. 1.2.1 - 4. 1.2.4 shows the abrupt shock jumps in all of the three plasma 
parameters - density, temperature and pressure — as the spacecraft crosses the bow shock. The 
spacecraft passes through the bow shock twice in one orbit interval: at both the inbound and 
outbound portions of each trajectories near the planet (Figure 4. 1.2.5). In both portions of the 
trajectory, all of the three solar wind plasma parameters generally remained steady at a certain 
value when the spacecraft was at a high altitude. The high altitude data suggests that the 
measurement has been taken outside the range where the solar wind-planetar/ ionosphere 
interaction takes place. As the spacecraft altitude decreases, we see the sudden jump in the 
plasma parameter values. We observed the bow shock formation from these jumps which 
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indicated the instant increase of each of the parameters. After the parameter values reach their 
maximum, the values gradually decrease back to, and occasionally becomes lower than, the 
initial values. This shows that the plasma contents beyond the bow shock have less direct 
contacts with the solar wind as the altitude continues to decrease toward the surface of the Venus. 
At such low altitude, relative to the bow shock location, many of the solar wind ions must have 
already deflected away toward the night-side of the Venus ionosphere or have become absorbed 
and neutralized in the upper atmosphere. 

Most of the plots from the 825 orbits showed a similar result, i.e. the plot showed the 
abrupt jump of the parameter values. The inbound plot showed such a jump after a steady initial 
high altitude value and the outbound showed the jump first and the steady high altitude value 
afterward (See Figure 4. 1.2.1). However, there were some abnormal cases as shown in Figures 
4. 1 .2.2 - 4. 1 .2.4. Figure 4. 1 .2.2 shows the unusual scattering of the parameter values after 
entering the bow shock. The shock shown in Figure 4. 1 .2.3 is so small that it is almost non- 
measurable. Furthermore, some of the shock jumps were unrealistically high as seen in Figure 
4. 1 .2.4. We suspect that the reason behind these abnormal cases are related to the considerable 
complexity of the solar wind-ionosphere interaction. Other important factor that must be 
considered include the direction and magnitude of the interplanetary magnetic field, the varying 
solar activity, and the direction of the ORPA instrument's gates for sample collecting relative to 
the solar wind flow direction. Also the position at which the spacecraft crosses the bow shock 
relative to the planet-sun direction influences the magnitude of the jump shock (Figure 4. 1.2.6). 
For example, if the measurements are taken at the center of the day-side hemisphere of the 
planet, the direction in which the interaction between the solar wind and the 
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Venus ionopause is almost perpendicular. Thus, the interaction there produces a stronger shock 
than that produced near the terminator. 


B-22 



Figure 4 . 1 . 2.1 OWPA Plot of How Shock Crossing Orbit 320 lnboutul 








Figure 4 . 1.2.2 OR PA Plot of Bow Shock Crossing Orbit 355 Inbound 
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Figure 4. 1.2.3 ORPA Plot of Bow Shock Crossing Orbit 365 Inbound 







Figure 4. 1.2.4 ORPA Plot of Bow Shock Crossing Orbit 334 Inbound 





Figure 4. 1.2. 5 The Orbit of PVO as viewed perpendicular to the orbit plane. 

Periapsis — solid dot, apoapsis — solid square. X’s are at 1 hour 
intervals before and after periapsis. The arrow denotes the direction of 
spacecraft 



The orbit of PVO as viewed from over the north pole of Venus to show 
the latitude covered by the near polar inclination of 74.4 degrees. 



OCNOTES OflfiK SIDE 


Source: Colin (1989) 
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Figure 4.1.2. 6 Diagram of the Solar Wind- Venus Interaction 
Shows: 

• Position of the Spacecraft Crossing of Bow Shock Relative to the Planet-Sun Direction 

• Magnetic Field Pile Up on the Day-Side Hemisphere of Venus 
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Source: Colin (1989) 
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4.2 M, pr vs Comparison 
4.2.1 Spreiter Model 

Spreiter and his coworkers have derived a model to describe the gas dynamics of the bow 
shock formation (1966). They developed a mathematical expression that relates the plasma and 
field parameters on either side of the bow shock. This expression has been derived under the 
assumption that the magnetic field induced by the solar wind itself was parallel to the flow, i.e., 
vH B. Although this parallel flow condition does not always hold to be true, this assumption 
simplifies the mathematical derivation allowing Spreiter and his coworkers to apply their gas 
dynamic model to the solar wind-planet interaction. 

Their findings were that the ratio of free streaming density, p„ to the density just behind 
the bow shock, p b , could be expressed as: 

p*/pjj = [ (y -1) M_ 2 + 2} / [( y + 1) M_ 2 ] ( y = ratio of specific heats; 

M_= upstream Mach number ) 

A similar one has been developed for the corresponding temperatures. They have developed 
these gas dynamic equations, initially , by applying them to the case of the solar wind-earth 
interaction. The solar wind interaction with the upper ionosphere of the earth is similar to the 
one for the solar wind-Venus interaction, in the extent that the terrestrial obstacle is spherical and 
the shock boundary is paraboloidal. However, as stated previously, the earth has a strong 
magnetic field that shields the terrestrial ionosphere/atmosphere from the solar wind. Combining 
this solar wind shield-off together with the assumption of parallel flow, they assume no entrance 
of the magnetic field into the picture. (Since the case of relative movement we have E = - v x B, 
if the solar wind velocity and the interplanetary magnetic field are parallel, then the E = - v x B 
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will be 0, allowing neglect of electromagnetic influence to be justifiable.) This absence of 
magnetic field influence partially justifies their use of only fluid equations, (i.e. use of only 
hydrodynamic equations instead of magnetohydrodynamic equations, which would be more 
appropriate to describe our system which involves the magnetic-field influence.) 

This model contains another simplification; namely, the assumption of a collision- 
dominated medium to be the interacting medium of the solar wind.. Unless modified 
appropriately, this assumption is not wholly applicable to our study which examines the nature of 
the solar wind, which is a collisionless plasma. This assumption on the nature of medium along 
with the assumption on neglect of no magnetic field force, the Spreiter model should be 
incomplete when applied to the Venus-solar wind interaction. In order for the Spreiter model to 
describe the solar wind interaction with a planetary obstacle like Venus more fully, a 
modification is needed that allows for the effects due to a collisionless plasma and mass loading. 
Nevertheless, the Spreiter model seems to work to a degree. Thus, the logical step for us to take 
then is to test the Spreiter model against the actual value of Mach number (the parameter 
involved in the Spreiter model) computed on the basis of the measured shock jump. This 
comparison should give us an idea of what scope of modification is necessary in order for the 
Spreiter model to become more complete. We suspect that their simplifying assumptions will 
become critical pivotal factors, when actual data has been compared to the values derived from 
using their model and we expected this to be the case for the PVO Venus data. Throughout this 
research, we will assume the atmosphere around Venus to be a monatomic gas. Hence, a value 
for the ratio of specific heats, y = 5/3 is used. 
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• Mach Number from the Spreiter Model: M^, r 

Mass density, p, in the Spreiter model is applicable specifically for a neutral proton- 

electron gas. Since our ORPA electron data is much more reliable and useful than the OPA 

proton data, we try to transform the Spreiter model into the form that substitutes the proton mass 

density with the electron number density. According to Tatraliyay and his coworkers (1992), we 

can safely make an assumption that the number of helium atoms counted from the PVO data 

would be negligible due to the overpowering proton abundance (the actual helium contribution 

might increase the density by 1-10%, which would not significantly influence the results.). The 

subscripts *> and b refer to upstream and downstream components, respectively, of the flow either 

side of the bow shock. The ratio of the mass density can be written as, 

pjp b = rL*m» / n b *m b = n_ / A> (since mass density, p = n*m and m_ = m b , 

the masses cancel out) 

resulting in, 

n. / n b = [ (y -1) M, 2 + 2] / [( y + 1) M_ 2 ] 

We now solve the above equation for M., 

n. ( y + 1) M. 2 - n b (y -l)M. 2 =2n b 
M. 2 [n.(y+1) - n b (y -l)] = 2n b 
M.= { 2n b / [n. ( y + 1) - n„ (y -1)] ] lp - 


We will express this Mach number as M, pr . 



• Mach Number from the PVO Data: 

The standard procedure to calculate the mach number is to use its exact definition: the 
ratio of velocity of the object to the speed of the sound, 

^def — ^jolar wind ^ Cjoaad* 

The speed of sound may be obtained as, 

c * = (Y P / P) 1/2 * [ Ynk(T e + Tp) / n(m e +m p ) ] 1/2 

( pressure: p = InkT = nk(T e + T p ) ) 

Again, using reasonable approximations, we want to transform the equation above in a way that 
the ORPA electron data is to be employed. 

Since m e « m p , m e + m p = m p . (1) 

Due to the proton temperature variation, we have to be careful when approximating the 
temperature values. Instead of taking a rough average, we have decided to retain two cases for 
our temperature values: T p maximum possible and minimum possible,. By doing so, we will be 
able to discuss the whole range of pressures at the solar wind. The electron temperature, T e , in 
the interplanetary medium typically ranges from 100,000 to 200,000 °K. The proton 
temperature, T p , on the other hand, ranges from about 20,000 °K to 200,000, i.e., 0.2 to 1.0 of Te. 
This information allows us to make a simplifying assumption: at T e max, T e = T p which leads to 
T e + T p = 2T e ; at T e min, T e + T p = 1 .2 T e . Therefore, our T e range used for this study becomes 
1.2 T e £ (T e + T p ) s 2 T t . (2) 

Combining conditions (1) and (2), and canceling n's, we have, 

Cj = [ yk(T e + T p ) / m p ] 1/2 = [ 2ykT e / m„ ] 1/2 (max) or [ 1.2ykT e / m p ] 1/2 (min) 
Putting this expression back into the Mach number equation, we get 
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M def = v, ( 2ykT e / m„ ) - ,/2 (max) or v,( laykT. / m, ) (min) x 
This Mach number is to be obtained from the ORPA data. We will express this Mach number as 

Mdef. 

4.2.2 Collecting Parameter Values for Comparison 

We need the following variable values to perform the comparison of and M^. 

— v„ T e with constants: (y, k, nip) 

M,,* — n,,, n_ with constants: (y) 

We took the random numbers of data from orbits 1-825. For each of these orbits, there could be 
two crossings of the bow shock by the spacecraft on the day-side hemisphere of the planet (see 
Figure 4. 1.2.5). First of all, T„ and n, are obtained from the ORPA electron data. Figure 
4.2.2. 1 shows the method we used in determining the parameter values from each orbit. T e and 
ru. are the average values of the horizontal steady line in the ORPA plot of temperature and 
density before the spacecraft crosses the bow shock. The density, n b , on the other hand, is the 
point where the approximate linear curve, fitted against the gradual decline in the parameter 
intersects the line parallel to the vertical axis that touches the edge of the shock jump. Also, it 
needs to be noted that we converted the unit of T e from electron volts in ORPA plot into degrees 
Kelvin to apply the values into the equation. 

T (°K) = [(1.60 x 10-' 9 Joules) / (1.38 x lO’ 23 Joules)] [°K/ev] T (ev). 

Second, the solar wind velocity data was obtained from the OPA data. (Although the 
proton density data from OPA is considered somewhat unreliable and temperature data not at all 
useful, its velocity data is safe.) As seen from the Figure 4.2.2.2, the solar wind velocity is a 
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rather stable parameter. In some orbits, they become almost a straight horizontal line with so 
small of a fluctuations in its 24 hrs. period, in which the spacecraft was set to complete for each 

orbital motion. Although the time resolution of the OPA data is still quite poor, this stability of 

the solar wind velocity allows us to use its average value as the velocity in the estimation of M^. 

4.2 3 Calculation and Comparison ofM^ r & M drf 

Next step is to enter all the values obtained into the equations below: 

M, pr = [ 2n b / (n„ ( y+ l) - n b (y -1)] ] l/2 
M def = v, [ (1.2 ~ 2) ykT e / m p ] ~ ,/2 

The Table 4.2.3. 1 contains the results from the computation of M spr and M def , for each randomly 
chosen orbit. We divided the M def columns into two which correspond with the two extreme 
cases of the electron- proton temperature relationships. Also, we calculated each Mach number 
for both inbound and outbound measurements. Although, the values from these two 
measurements from one orbit differed slightly, their result came out to be quite similar. Thus, we 
only plotted the M^-M^ comparison graph only for the inbound, assuming the result from the 
outbound would be similar. When we see the Figure 4.2.3. 1 which shows the M spr plotted along 
with M def , we noticed, that except the two eases out of the thirty cases (orbit 40 and 90) we 

calculated, M spr was smaller than M^. Beside the two abnormal points,, the relationship, M spr < 

M det -, was consistent. When we take the plasma parameter values under the condition of T e = T p , 
the average value of was 2.71 with its maximum 4.27 (excluding the two abnormal cases) 
and its minimum 1.14. Similarly, the average value of M def , on the other hand, was 6.93 and its 
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Tabic 4.2.3. 1 Computed Results of M spr and M 
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minim u m 2.28 and its ma ximum 1 1.06. We noticed that does not bounce around as much as 
M^. This could be due to our lack of knowledge of exact proton temperature, since we 
substituted the electron temperature which shows a considerable variation. 

As noted above, the Mach number derived for orbit 40 and 90 showed abnormal behavior 
compared to the rest of the Spreiter Mach number data: M^ (#40) = 7.97 and M^ (#90) = 7.55. 
As we studied the original ORPA plots of plasma parameter versus the spacecraft altitude, we 
found that the orbit #90 had a relatively low electron temperature with a small jump (see Figure 
4.2.3 .2) and that the orbit #40 had an extremely high density jump (Figure 4.2.3 .3). These 
abnormal characteristics the original plasma parameter showed in both orbits might have led to 
the contrasting result to the general M spr < M det - relationship. We also noticed that the value of the 
M def becomes closer to M spr in these abnormal situations, when we assume the proton temperature 
to be considerably smaller than the electron temperature, i.e., T e » T p , instead of T s = T p 
condition. For example, in the case of the orbit # 40, M def (T e » T p ) = 6.26 whereas M def (T e - 
T p ) = 4.43 widens its gap against the M det - = 7.97. Again, our uncertainty in the proton 
temperature may be effecting the result in these abnormal situations. It is also possible that the 
magnitude and the direction of the interplanetary magnetic field produced an abnormally strong 
type of interaction. 

4.2.4 The Possible Explanations to the M, pr -M def Variance 

After carefully examining the result of M spr -M def comparison, we believe that we have 
verified that the Spreiter model is incomplete as far as its predictability of the bow shock 
formation for planets similar to Venus. M spr did not agree well with M deI -, -- M 3pr was generally 
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Figure 4.2.3.2 ORPA Plot of Bow Shock Crossing Orbit 90 Inbound 
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Figure 4.23.3 ORPA Plot of Bow Shock Crossing Orbit 40 Inbound 
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found to be smaller than by the average ratio of about 1 to 2.6, respectively, which was a 
big difference, even though they both were describing the same phenomenon. Thus, we then 
searched for the possible explanations to why the referred Mach number derived using the 
Spreiter model differed so much from the actual Mach number computed on the basis of its 
definition. 

We suspected the most crucial vulnerability of the Spreiter model to lie in its initial 
assumption that neglects its magnetic field influence on the solar wind plasma. The base of the 
Spreiter model came from its application on earth, which has a definitely different magnetic field 
than that of Venus. Let us discuss the differences in the nature of the solar wind interaction 
with both planets, earth and Venus, and simultaneously refer back to the Spreiter model 
feasibility. By doing so, we may see what kinds of modification the Spreiter model needs in 
order for it to be applicable to other planetary- obstacles with different planetary magnetic fields. 

The terminology, “hard” and “soft” are sometimes used in the field of plasma physics to 
classify the solar wind obstacles depending on how rigidly each planet can shield off the solar 
wind entrance into its ionosphere region. The Spreiter model has been created under the 
restriction that the solar wind is interacting with a "hard" obstacle, such as the planetary magnetic 
field of the earth. That is, the earth is considered to be a "hard" obstacle because of its strong 
magnetic field. If the solar wind is able to interact with the ionosphere and upper atmosphere of 
a planet, the obstacle presented by the planet will likely be “soft” because of mass loading 
(discussed previously) and other factors. As mentioned before, Venus' magnetic field is almost 
negligible or undetectable, thus, Venus is considered to be a "soft" obstacle by allowing the solar 
wind to interact directly with its atmosphere in its ionopause region. The Venus' bow shock and 



magnetosheath are located much closer to the planet than that of the earth. The word, “soft” 
seems appropriate, since the Venus is like a "soft cushion" that is unable to resist the external 
impact on its surface. On the other hand, the earth and its strong magnetic field appear hard like 
"concrete," guarding off the impact of the solar wind, allowing very small influence by the 
intruder, the solar wind, in its ionosphere. 

When the “soft” obstacle lets the solar wind into its ionosphere, various 
magnetohydrodynamics reaction may take place. As mentioned earlier, one such interaction 
involves the charge exchange between the solar wind protons and the atmospheric neutral oxygen 
ions. The charge exchange results in the pick up of these ions by the solar wind electric field 
(E = -v x B) ala E x B drift and eventually to a mass loading effect. The slowing down caused 
by the mass loading of upstream plasma near the bow shock region causes formation of the 
shock jump, which is usually smaller than that due to a strong planetary magnetic field. 

It seems reasonable to assert, therefore, that the Spreiter model becomes inapplicable 
when the obstacle is a “soft” obstacle. The more dominant the magnetic field influence is 
allowed to become, due to the inability of the planet’s defense against such intervention, the 
more complex its influences become on the plasma interaction around the planetary obstacle. 
Since the Spreiter model does not take the magnetic field influence into account, the softer the 
obstacle becomes, the more the model appears to show its vulnerability of not being able to 
respond to all of the MHD reactions taking place upstream and in the vicinity of the bow shock. 
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43 Pick-Up Ion Fraction 

43.1 Bow Shock Comet Model for Ion Pick-Up 

One of the notable MHD reaction that may be taking place in the solar wind-Venus 
interaction is that of a mass loading effect. Mass loading effect is less likely to occur on the 
“hard” obstacle such as earth because the planetary magnetic field prevents the solar wind from 
getting close enough to the atmosphere/ionosphere. As noted previously, the discrepancies we 
have observed when comparing the plasma parameters, in the forms of M spf and M def , may be 
due to the Spreiter model’s inability to include the MHD effect such as mass loading in its 
expression. However, we are not sure to what extent this mass loading effect influences or 
causes the bow shock formation. Hence, the next logical step is to try to evaluate the 
significance of this mass loading effect. In order to measure the significance of such an effect, 
we may estimate the amount of ions that get picked up as the solar wind interacts with the 
planet’s ionosphere. The amount of ion pick up occurring should be directly correlated to the 
relational impact of the mass loading on the bow shock formation. 

There is a model developed by Cloutier (1982) which describes the fractional ratio of ion 
pick up taking place during the solar wind-comet interaction. This model has been developed 
as part of the study of cometary bow shocks. A comet is another obstacle that interacts with the 
solar wind. Although there are similarities, the comet also differs from the earth and Venus in 
various wavs. A comet has an extremely small nucleus, nevertheless, at 1 A.U. of the sun the 
virtual absence of gravity allows its atmosphere to be very large, much larger than that of earth or 


B-44 



Venus. Because a comet does not have a magnetic field, but its atmosphere is so large, it is 
considered to be “extremely soft” Venus is a similar obstacle to a comet in the context that it is 
incapable of shielding off the solar wind from the ionosphere, but because of its strong gravity, 
the shock interaction is much closer to the planet The comet’s bow shock formation is almost 
exclusively due to mass loading pick up, whereas that at Venus must also be due in part to 
magnetic field pile up. Furthermore, the use of Mach number in the equation used in Cloutier’s 
model is another aspect that makes this model a good candidate to be applied in our study. 

Hence, it seems logical to consider the Cloutier model as another extreme to apply to the PVO 
data to help better understand some of the complexities of the Venus bow shock formation. 

Before introducing the model, we need to understand the movement of the ion when 
created upstream of the nucleus of a comet. As similar to the case of Venus, the mass loading 
occurs near the bow shock. This mass loading takes place due to the ions created upstream in the 
flow of the solar wind by whatever process that involves the ionization. The ion finds itself in 
the electric field in the instant it becomes an ion, which is expressed as E = -v x B accordingly 
with the Faraday effect for the geometry below: 



Then, the particle starts its drifting motion (as shown in Figure 4.3. 1.1) in the direction of the 
flow with a velocity given by v D . The drifting motion of a particle under magnetic field 
influence is expressed as 
v D = E x B / B : 

B -45 



Figure 4.3. 1.1 Drift Motion of Pick Up Ions Due to the Solar Wincl-Venus Interaction 



Source: Colin (1989) 


which ends up with the solar wind velocity, v, 

v B = ExB/B 2 = ExB/B 2 = -(vxB)xB/B : 

= B x (v x B) / B 2 = v B • B / B 2 - B (v *B) / B 2 
= v 

The fact that, v D = v, shows the conservation laws of motion being held within the plasma 

during the interaction between the solar wind and the comet, which will be used in Cloutier’s 

model discussed later. Note that when 0 + is formed, it has a velocity of only several kilometers 
per second The action of E x B / B 2 to bring 0" up to solar wind velocities must slow down the 
solar wind. 

Wallis (1973) has shown that the bow shock formation is linked to the condition in which 
the modified mass flux p V satisfies the equation: 
p'v’ = P , v, ( y 2 If -1) 

where y is the ratio of specific heats (used also in Spreiter model), p s is the upstream 
undisturbed solar wind mass density, and v, is the solar wind velocity. Wallis transforms this 
model into a form in which the ion pick up fraction is calculatable. The modified mass flux, 
p’v’, is the resultant flux after the additional flux, p v, has been added due to the created ions 
during the mass loading process as shown in the diagram below. 

— — > --> 

P , v, 11“ P* V 

S = p v added 

However, he oversimplifies the phenomenon by only applying the conservation of momentum 
flux. Years later, another plasma physicist, Cloutier made improvements on the Wallis’ model 
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by adding the mass flux and the energy flux into the picture (the Wallis* model had only included 
the momentum flux). Thus, the number of the needed flux equations become three and they are 
expressed as: 

Mass p jVf +S =pV (3) 

Momentum p’ v’ 2 + p’ = p , v, 2 + p , (4) 

Energy [ 1/2 p V 2 + (y / ( y -1)) p’] v’ = [ 1/2 p’ v’ 2 + (y / ( y -1» pj v* (5) 

The variable q is employed to simplify the process of solving the equations and also to perform 
the derivation in terms of the ratio of the velocities, 

T[ = V, / v\ 

Combining the equation (3), (4) and (5), we obtain the equation for q, 

T1 = 1/2 {(l+q 0 ) +/- [(1- q 0 ) 2 - 4 or q 0 ] "*} 

where 

m [( y + 1) M 2 ] / ( 2 + ( y - 1)J 
M 2 = I v, / ( y p, / pj 1/2 ] 2 = p, v, 2 / y p 5 

o =S/p,v, 

Here, we face the possibility of q becoming an imaginary number. The imaginary number in the 
physical expression indicates that the system is in an unstable condition, hence, the system is 
likely to go through a non-smooth transition. Thus, if tj is imaginary, we have a good possibility 
of having a jump in the parameters, which leads to the formation of the shock. The satisfaction 
of the condition, 4 cr q 0 s (1- q 0 ) : , is necessary for the shock, i.e., a shock forms if 
a > (1- qj 2 / 4 ru. (6) 
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The right-hand side of (6) will turn out to be, according to the Cloutier’s model, the threshold 

point for the shock formation. Hence, we call it 
= (1- r\ J 2 / 4rj 0 . 

This model enables us to calculate the pick up ion fraction, a, which is related to the ratio 
of mass density p’/p, obtained through the manipulation of the equations introduced above. 

4.3.2 Calculation of Pick-Up Ion Fraction 

Although the typical Mach number, M = 2, for the solar wind-comet interaction is 
usually used for the calculation of ion pick up fraction for the Cloutier model, we will use the 
average Mach number we obtained using the Spreiter model, to estimate the ion pick up fraction 
in the case of Venus. We will first take the average M spr of 2.71 with the ratio of specific heats of 
5/3. 

Thus, for M = 2.71 and y = 5/3, 

T] 0 =[(y+l)M 2 ]/[2 + (y-I)l = 2.8399 

thus = (1-t 1o ) 2 /4t 1o = 0.2980 

and q = v, / v* = 1/2 {(1+q.) +/- [(1- qj 2 - 4 a r| 0 ] 1/2 } = 1.9228. (7) 

Using a = S / p, v, to replace flux S, 

p' V = PjVj +S => p’ v’ = PjVj + apjV, = p.Vjfl+a) 

Hence, 

p' v’ = p, v, (1 + 0.2980) = 1.2980 p, v, (i.e., 0.2980 p s v s added) 

p'/p, =1.2980 v’/v. 
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= (1.2980) (1.9228) 
= 2.4956 


( since t| = v,/ v* from (7) ) 


( 8 ) 

Putting the mass density equation in terms of the ion pick up fraction value, a, (0 < a < 1), 

p' = n (16 a + (1- a) ] m p (9) (16mp for mass of oxygen) 

P, =am, (10) 

By equating the ratio of mass density before and after the flux addition, we obtain an equal form 
of mass density ratio as (8), hence, 

p' /p, = In {(16 a + (1- a) } m p ] / n m p 

= 16 a + (1- a) 

= 2.4956 

From the relationship above we compute a, 

16 a + (1- a) = 2.4956 

a = 0.09971 

Finally, we obtained the value for the pick up fraction of 16 mp ions (oxygen ions) using the 
average Mach number from the Spreiter model. 

a ave = 9.97 % 

Similarly, we computed two other ion pick up fraction values to determine the range of a by 

taking the minimum and maximum M, pr values resulting in values, respectively of 
a = 0.798 % and 23.63 %. 

From the fractional values obtained above, we can discuss the significance of the pick up 
of oxygen ion (O') taking place near Venus ionosphere. The occurring fraction of 9.97 % is 
relatively high. Almost one out of ten solar wind ions coming in will give its charge away to the 
neutral oxygen atoms that flowed outwardly from the Venus atmosphere. Hence, we argue that 
the ion pick up is a factor significant enough to be included into the bow shock formation 
models. Although applying the comet model to the Venus problem is oversimplifying the 
physics, we suggest that this result is meaningful. It provided us the idea about the importance of 
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MHD influences in the bow shock forming medium, particularly of the ion pick up which is the 
major cause of the mass loading of the solar wind. 

We also noticed the wide range the fractional value, a, lies in, i.e., 

<*.«/«»»= ( 23.63 %) / ( 0.798 %) = 29.6 = 30 

While (M,,, max) / (M^ min) = 3.6). This shows that the a is very sensitive to the change in the 
Mach number. At the same time, we know that the inferred a must be strongly connected to the 
nature of the interplanetary magnetic field embedded in the solar wind, since its piling up at the 
ionopause could modify the amount of charge exchange that takes place. It seems logical, 
therefore, to expect a relationship between the Mach number and the strength and orientation of 
magnetic field, and also that such a relationship would vary as the a value fluctuates. We 
suspect that the high variation in the rate of ion pick up occurring is likely due to the magnetic 
field variation and also to the changing solar activities. 

Two conditions, one from the magnetic field variation and the other from the changing 
solar activities, seemed to possibly influence the resulting ion pick up fraction. Let us discuss 
the influence from these two factors: the direction of magnetic field and the solar activity cycle. 

If the interplanetary magnetic field is perpendicular to the path of the solar wind, the 
highly conducting state of the Venus ionosphere could lead to a pile up of the interplanetary 
magnetic field on the day-side hemisphere of the planet.(see Figure 4. 1.2. 6) The magnetic field 
pile up could lead to the creation of a “hard” obstacle out of a “soft” obstacle, since the piled up 
magnetic field can behave as a shield against the solar wind intervention into the ionosphere of 
Venus. If the magnetic field achieves such a condition, the solar wind must deflect away from 
the planetary obstacle, and mass loading due to ion pick up is less likely to occur under such 
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conditions. 


The other case is solar cycle influence on the ion pick up. When the sun is active and 
producing intense UV radiations, this could cause more ionization of oxygen atoms upstream of 
the planet’s ionosphere/atmosphere, and increase the amount of mass loading. (See Figure 
4.3.2. 1 which shows the solar cycle effect on Venus' bow shock.) However, at such times the 
interplanetary magnetic field may be stronger. It is not clear which of these two competing 
factors are more important 
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Figure 43.2. 1 Solar Cycle Effects on Venus’ Bow Shock 




Year 


Source: Colin (1989) 



. Conclusion and Summary 


. By utilizing plots of the bow shock jump using the ORPA electron data, we have tried 
to test the feasibility of Spreiter’s gas dynamic model which attempts to describe the nature of 
the bow shock formation. The comparison between the inferred mach number using the Spreiter 
model and the experimental data from its definition showed the relationship of < M def to be 
true most of the time. The Mach number calculated using the Spreiter model had its average of 
2.71 and ranged from 1.14 to 4.27 (excluding two abnormal events), whereas the Mach number 
from its definition had its average of 6.93 and ranged from 2.28 to 1 1.06 . 

We argue that the reason behind such a disagreement between the two values is the 
oversimplification on the initial assumption the Spreiter model makes. The Spreiter model was 
initially developed for a “hard” obstacle such as the earth. The earth has a strong magnetic field 
that can shield the solar wind from the planet’s ionosphere/atmosphere. However, for the case 
of Venus, there is no planetary magnetic field, and, therefore, the solar wind freely enters into its 
ionosphere and starts to interact with the neutral atoms flowing upstream of the planet. This solar 
wind’s entrance into outer atmosphere of the planet creates various complexities when one 
attempts to predict the bow shock formation, which the Spreiter model is not able to do in the 
case of Venus. 

One of the complexities the solar wind- Venus interaction faces and Spreiter fails to 
include in his model is that of the mass loading phenomenon. Mass loading at the solar wind- 
planet interaction points slows down the plasma flow in the region. We tried to estimate the 
impact this has on the plasma parameter fluctuation and, consequently, to the bow shock 
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formation. We used Cloutier’s model developed for the solar wind-comet interaction which 
include the expressions which allows us to calculate the ion pickup fraction, which is directly 
related to the mass loading. We reasoned that if the ion pick up fraction is great, then mass 
loading factor is significant enough to be included into the basic assumption which the predictive 
model on bow shock formation would be based on. 

The results from such a computation using the Cloutier’s model was that the fraction of 
16 nip ions (oxygen ions), that were picked up had its average of 9.95% and ranged from 0.78% 
to 23.63%. We argue that these fractional values are significant enough, and therefore mass 
loading effects should be considered when constructing a model for bow shock formation, 
particularly, for the planetary obstacles that are “soft,” i.e., that are vulnerable to the solar wind 
intervention into the planet’s ionosphere region. However, this result takes only charge 
exchange and mass loading into account, neglecting the effect of interplanetary magnetic field 
pile up in the day-side ionosphere of Venus. This is clearly needed in order to more accurately 
estimate the amount of charge exchange production and subsequent pick up of Q-K 

Clearly, there are various complex MHD actions taking place in the bow shock forming 
region. We showed our concerns on the neglect of these significant factors in the already- 
existing bow shock models. However, we also faced the difficulty in modifying these plasma 
dynamics models due to its overwhelming complexity and were forced to neglect important 
factors such as the direction of magnetic field and solar activities. However, this research 
confirms the need for a continuation of careful investigations on this uniquely complex solar 
wind- Venus interaction, possibly employing the OPA-ORPA data available from the PVO 
mission which still may offer us the various uncovered truths about the complex nature of 
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plasma. Moreover, this study also should draw the attention of future researches in that a model 
with the various MHD factors involved needs to be developed in order for us to further our 
understanding on the curious phenomenon of bow shock formation and the complex nature of 
plasma dynamics occurring at Venus. 

We will close with the following quote from Ten Years of Discovery, the publication 
completed by the members of PVO mission team, since it seemed to echo with this Venus bow 
shock study and many other scientific efforts all over the world trying to understand the 
complexity of the universe. 

. .Pioneer Venus discoveries have to a large degree answered the simple questions 
about the Venus ionosphere that were posed before encounter. But, in the process of 
answering he simple questions, the PVO investigators have become aware that Venus is 
much more complicated than our early questions implied. Thus 1 0 years later we sill 
have many questions, and the answers to our original questions have been incomplete in 
many cases. They were the right initial questions, of course, but nature is always more 
complex than one imagines from a position of ignorance. Our direct contact with this 
new world has brought with it many surprises and a new set of question: questions 
concerned not so much with what the Venus ionosphere is like, but what causes it to be 
that way. . . 
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